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' Abstract. How does the submm galaxy population detected by SCUBA fit 

I into galaxy evolution? How do these rare starbursting systems, which contribute 

significantly to high redshift star-formation, relate to other galaxy populations? 
Deep radio observations have been most useful for studying these systems, but 
still leave a significant fraction of the population unidentified. Now with the 
deep IRAC and MIPS images from the GOODS Spitzer Legacy program and a 
re-analysis of the deep radio data, we are able to identify counterparts for a large 
fraction of SCUBA galaxies in GOODS-N. All of these counterparts are detected 
^ ■ by Spitzer . Given the vast multi-wavclcngth data in this field, we can study the 

' spectral energy distributions (SEDs) of these systems and determine what is fu- 

, eling their intense infrared luminosities. A rest-frame composite optical-to-radio 

'sj" ■ SED for all spectroscopically identified submm sources shows that the average 

' SCUBA galaxy is consistent with models of ultraluminous starburst galaxies, 

, although cooler than those observed locally. Because of this, the submm flux 

alone consistently overestimates Ljr when using spectral templates which obey 
' the local ULIRG temperature- luminosity relation. The wide range of 24/850 

, micron flux ratios as a function of redshift indicates the presence of strong mid- 

Q^' IR features, to be confirmed with deep IRS spectroscopy. The IRAC colours of 

the submm systems provide useful redshift constraints, since, at these redshifts, 
IRAC samples the stellar bump. The Spitzer photometry of this large sample of 
submm galaxies has allowed us to put constraints on many of the outstanding 
issues in submm astronomy. 
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■ 1. Introduction 



The Submillimetre Common User Bolometer Array (SCUBA, Holland et al. 1999) 
on the James Clerk Maxwell Telescope (JCMT) is efficient at finding high red- 
shift dusty galaxies in blank-field extragalactic surveys. These galaxies are Ul- 
traluminous, and are known to contribute significantly to the star-formation 
budget at high redshift (e.g. Lilly et al. 1999; Chapman et al. 2005). However a 
general understanding of the submm population is severely limited by the small 
number of known SCUBA sources, the large beam size at 850 fj,m, and their 
faintness at other wavelengths. Enormous progress has been made using radio 
observations to localize the submm sources (e.g. Ivison et al. 2000; Small et 
al. 2000; Barger et al. 2000). At least 1/3 of currently known 850 /xm galaxies 
have no radio counterpart in some of the deepest radio observations and there- 
fore deep observations at other wavelengths are crucial. No other SCUBA field 
has the extensive multiwavelength coverage of the Great Observatories Origins 
Deep Survey (GOODS) Northern field. In addition to the VLA radio imaging 
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Figure 1. Multi- wavelength images of four submm counterparts, from left to 
right: ACS zsso; IRAC 3.6 //m; MIPS 24 /xm; and radio 1.4 GHz. All images 
are 20 x 20 arcsec and the black circle is centred on the SCUBA position 
with a 7 arcsec radius. The smaller black/white circle indicates the submm 
counterpart. Submm galaxies are quite faint at other wavelengths, but deep 
MIPS data are extremely valuable for identification. 



(Richards 2000; Morrison et al. in prep), deep Hubble Space Telescope {HST) 
imaging in four bands (Giavalisco et al. 2004) and the Chandra 2 Msec image 
(Alexander et al. 2003), GOODS boasts the deepest observations taken in the 
infrared with IRAC and MIPS on the Spitzer Space Telescope (Dickinson et 
al. in preparation). GOODS-N is also home to a large amount of SCUBA data 
taken by several groups (Hughes et al. 1998; Barger, Cowie, & Richards 2000; 
Borys et al. 2003), which have been combined into one large 850 ^m. 'super-map' 
(Borys et al. 2003, 2004; Pope et al. 2005). 



2. Identifications 



We have identified 35 secure 850 /xm sources in the GOODS-N submm super-map 
(Borys et al. 2003; Pope et al. 2005). Using all multi- wavelength data, unique 
counterparts are found for 94 per cent (33 out of 35) of these submm sources. 
26 out of 35 are detected above 2>a in the new Morrison et al. reduction of the 
1.4 GHz data and all of the counterparts are detected with IRAC and MIPS. 
Fig. n shows multi- wavelength images of four of our submm sources. 



Using Spitzer to study submm galaxies 



3 




Figure 2. Left: Cumulative distribution of 24 /xm flux for submm counter- 
parts in GOODS-N. Right: The S'24/5'850 ratio as a function of redshift for 
submm sources in GOODS-N. The dashed curve is an SED template from 
Chary & Elbaz 2001 (CEOl) models, the solid curve is a modified CEOl tem- 
plate which includes additional extinction from Draine (2003), the dash-dot 
curve is the observed SED of Arp220 and the dash-dot-dot curve is Mrk231, an 
IR-luminous AGN. The shaded region represents the region where the 9.7 /im 
silicate feature passes through the 24 /zm passband. 



3. MIPS 24 fim properties 

Both SCUBA and Spitzer are able to select distant star-forming galaxies which 
contribute significantly to the cosmic infrared background. In blank field sur- 
veys, SCUBA is primarily sensitive to the thermal dust emission from ultra- 
luminous galaxies at high redshift. While Spitzer also detects these galaxies, its 
sensitivity is such that it detects galaxies over a wider range of luminosity and 
galaxy type. The left panel of Fig. |21 shows the MIPS flux distribution for our 
850 /im selected sample. The vertical dashed lines indicate the 3a limits of sev- 
eral other Spiteer surveys: GTO shallow survey of the Bootes field (E. Le Floc'h 
priv. communication), FLS, SWIRE (D. Farrah priv. communication) and the 
GTO deep survey of the Lockman Hole East region (Egami et al. 2004). The 
GTO shallow survey of the Bootes field is deep enough to detect 25 per cent of 
the SCUBA galaxies. With the MIPS depth of the First Look Survey (FLS^), 
or SWIRE (Lonsdale et al. 2004), we would detect 40 per cent of our submm 
counterparts. 

In the right panel of Fig. |2 the S'24/5'850 ratio is plotted as a function of 
redshift for the sample of submm sources in GOODS-N as compared to several 
models. The generally lower S24/S850 ratios for the submm sources suggest that 
they have higher levels of extinction, especially in the mid-IR, than the local 
ULIRGs which were been used to construct these models. In particular, the 
silicate absorption feature at 9.7//m may be attenuating the 24 //m flux. 



^ http://ssc.spitzer.caltech.edu/fls/extragal/spitzer. html 
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Figure 3. Left: IRAC-based colour-magnitude diagram for GOODS-N. 
Larger symbols are our submm sources, while the small black dots are field 
galaxies. The squares are at low redshift (z < 1.5), the circles are at high red- 
shift (z > 1.5) and the triangles have unknown redshift. Crosses are sources 
which contain an AGN, as indicated by the presence of hard X-rays. Right: 
Composite SED for the GOODS-N submm sources with spectroscopic red- 
shifts. IRAC through radio photometry data are plotted. The solid curve is 
a modified CEOl model for a cool starburst galaxy, where we have applied 
additional extinction using the Draine (2003) models to minimize scatter rel- 
ative to the data points. The dash-dot curve is a CEOl model with the same 
luminosity but a different peak wavelength. This template does not have any 
additional extinction and is a poor fit to the mid-IR and submm data points. 



4. Spectral energy distribution 

At the redshifts of these submm galaxies, the IRAC channels are sampling the 
rest-frame near-IR, which is sensitive to direct emission from stars. If these 
galaxies are not dominated by an AGN at these wavelengths then the IRAC 
colours should show an indication of the 1.6 /xm local peak of the fi, stellar 
SED. In the left panel of Fig. |21 we plot an IRAC-based colour-magnitide di- 
agram, specifically S^s/Sss as a function of Ss.s, for the submm sources and 
the GOODS-N field galaxies. The redshift of the submm source is indicated by 
the symbol type. We see that z ~ 1.5 provides a clear transition in this IRAC 
colour for the submm sources, due to the filters moving over the stellar peak. 
The sources which lack redshift information are all consistent with the colours 
of the higher redshift submm sources. There seems to be no strong colour dis- 
tinction between the sources with and without a hard X-ray detection, although 
the former are generally brighter at 5.8 ^m. 

The right panel of Fig. |31 shows a composite SED for the submm sources in 
our sample which have spectroscopic redshifts, as compared to several models. 
The main conclusion to be drawn from this plot is that high redshift submm 
galaxies appear to have different SED shapes from those observed locally for 
galaxies of similar luminosities. A major contribution to this difference could be 
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Figure 4. Infrared luminosities and SFRs for submm sources. The bottom 
axis show the luminosities estimated using only the 850 fiia flux and the red- 
shift. The left-hand axis is the estimated IR luminosity using data at 24 /im, 
850 /im and 1.4 GHz. 



the temperature of the dust, but it is important to realize that this just means 
an SED peaking at longer wavelengths. The best-fit submm SED (solid curve in 
Fig. 131) peaks at 100 //m (T ~ 30 K), while a typical local ULIRG template from 
CEOl (dash-dot curve in right panel of Fig. ^ of the same luminosity peaks at ~ 
85 fxm (T ~ 34 K). This result emphasizes the strong selection effects, both locally 
and at high redshift. Our knowledge of local ULIRGs is not immune to selection 
effects, since it is dominated by results from IRAS which may preferentially select 
galaxies with warmer SEDs (peaking at shorter wavelengths). Similarly at high 
redshift, 850 /xm imaging preferentially selects galaxies which peak at longer 
wavelengths, and therefore comparing the two samples becomes problematic. 



5. Infrared luminosities 

Using the CEOl models, we have estimated the IR luminosity by fitting to the 
templates: 1) the 850 ^m flux and the redshift, and 2) the 24 ^um, 850 /xm and 
1.4 GHz flux and the redshift. Fig. compares the infrared luminosity (LIR) 
as estimated these two ways, by integrating from 8-1000 iJ^m under the best- 
fit modified CEOl template. Filled circles and squares denote submm sources 
with spectroscopic and photometric redshifts, respectively. The crosses denote 
sources which contain an AGN, as indicated by the presence of hard X-rays. The 
top axis shows the corresponding star formation rate (SFR) using the expres- 
sion in Kennicutt (1998). Using templates based on local ULIRGs with a fixed 
luminosity-temperature relation (CEOl), the submm flux will typically overes- 
timate the luminosity. This effect is most dramatic for less luminous submm 
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objects which tend to he at lower redshifts where the 850 /xm window samples 
even further from the IR SED peak. 

The median IR luminosity (using radio, submm and infrared flux) for our 
sample is Lg-iooo/nm = 6.0x IO^^Lq, which corresponds to a SFR of 1100 MQyr"^ . 

6. Summary 

Using the deep Spitzer Legacy images and a re-rcduction of the VLA radio data, 
we find secure counterparts for a large sample of submm galaxies in GOODS-N. 
The 24 /xm to 850 /nm flux density ratio as a function of redshift is low compared 
to models, suggesting that the silicate absorption feature at 9.7 /um may be at- 
tenuating the 24 iim flux. A composite rest-frame SED shows that the submm 
sources peak at longer wavelengths than those of local ULIRGs. This demon- 
strates the strong selection effects, both locally and at high redshift, which may 
lead to an incomplete census of the ULIRG population. Wc determine the IR 
luminosity Lir(8 — lOOO/xm) by fitting models and find that the submm fiux 
alone consistently overestimates Lir when using spectral templates which obey 
the local ULIRG temperature-luminosity relation. The shorter Spitzer wave- 
lengths sample the stellar bump at the redshifts of the submm sources, and we 
find that the Spitzer photometry alone provides a model independent estimate 
of the redshift. More details on this work will be presented in Pope et al. (2006). 
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